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We have theoretically investigated macroscopic quantum tunneling (MQT) and the influence of 
nodal quasiparticles and zero energy bound states (ZES) on MQT in s-wave/ d-wave hybrid Joseph- 
son junctions. In contrast to d-wave/d-wave junctions, the low-energy quasiparticle dissipation re- 
sulting from nodal quasiparticles and ZESs is suppressed due to a quasiparticle-tunneling blockade 
effect in an isotropic s-wave superconductor. Therefore, the inherent dissipation in these junctions 
is found to be weak. We have also investigated MQT in a realistic s-wave/d-wave (Nb/Au/YBCO) 
junction in which Ohmic dissipation in a shunt resistance is stronger than the inherent dissipation 
and find that MQT is observable within the current experimental technology. This result suggests 
high potential of s-wave/d-wave hybrid junctions for applications in quantum information devices. 
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I. INTRODUCTION 

Since the experimental observations of macroscopic 
quantum tunneling (MQT) in YBCO (YBaCuO) grain 
boundary^ and BSCCO (BiSrCaCuO) intrinsioM^ 
Josephson junctions, the macroscopic quantum dynam- 
ics of high-T c (i-wave junctions has become a hot topic 
in the field of superconductor quantum electronics and 
quantum computation. Recently, the effect of low-energy 
quasiparticles, e.g., the nodal-quasiparticles and the zero 
energy bound states (ZESs)&£ on MQT in <i-wave junc- 
tions have been theoretically investigated] 10 i 11 ' 12 i 13 i 14 It 
was found that, in c-axis type junctions^ the suppres- 
sion of MQT resulting from the nodal-quasiparticles is 
very weak i 10 ' 11 This result is consistent with recent ex- 
perimental observations In the case of in-plane 
type d-wave junctions^ however, the ZESs give a strong 
dissipative effecti 12 : 13 ' 14 Therefore, it is important to 
avoid the formation of ZESs in order to observe MQT 
with a high quantum-to-classical crossover temperature 

On the other hand, recently, quiet qubits consisting of 
a superconducting loop with an s-wave/d-wave (s/d) hy- 
brid junction have been proposedi 17 ' 18 : 19 In quiet qubits, 
a quantum two level system is spontaneously generated 
and therefore it is expected to be robust to decoherence 
from fluctuations of the external magnetic field i 17 ' 18 : 20 
However, the influence of the low-energy quasiparticle 
dissipation due to nodal-quasiparticles and ZESs on s/d 
quiet qubits is not yet understood. Therefore, it is impor- 
tant to investigate such intrinsic dissipation effects on the 
macroscopic quantum dynamics in order to realize quiet 
qubits. 

In this paper, motivated by the above studies and the 



recent phase-sensitive spectroscopy experiments in s/d 
junctions ; 21 : 22 we will discuss the application of a generic 
MQT theor y 23 : 24 to s/d hybrid Josephson junctions [see 
Fig. 1(a)]. In contrast with d/d junctions , 10 i 12 : 13 we will 
show that the influence of the low-energy quasiparticle 
dissipation is suppressed by virtue of the small but finite 
isotropic gap in the s-wave superconductor. Therefore, 
weak quasiparticle dissipation is anticipated in such junc- 
tions. We will also discuss an extrinsic Ohmic dissipative 
effect in a real s/d (Nb/Au/YBCO) junction 21 .^ and 
show that high T* , comparable to high-quality s/s junc- 
tions, is expected. These results clearly indicate the ad- 
vantage of s/d junction for qubit applications with longer 
coherence time. 



II. THEORY OF MACROSCOPIC QUANTUM 
TUNNELING 

A. Macroscopic quantum tunneling in 
s-wave/d-wave junctions 

In the following we will calculate the MQT rate in s jd 
junctions with a clean insulating barrier and without ex- 
trinsic dissipation, e.g., an Ohmic dissipation in a shunt 
resistance. Extrinsic dissipative effects in an actual junc- 
tion will be discussed later. The partition function of a 
junction can be described by a functional integral over 
the macroscopic variable (the phase difference <^) ) 24 i 25 : 26 
i.e., 

Z = /^)exp(_^l). (1) 
In the high barrier limit, i.e., zq = mw^/t^kp 3> f (m is 
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FIG. 1: (Color online) (a) Schematic of the in-plane s-wave/d- 
wave hybrid Josephson junction, a and b denote the crys- 
talline axes of the d-wave superconductor, and \ is the mis- 
match angle between the normal of the insulating barrier (I) 
and the a-axis. (b) Potential U((f>) v.s. the phase difference 
<j> between two superconductors. uo p is the Josephson plasma 
frequency of the junction. 



the electron mass, Wo is the height of the delta function- 
type insulating-barrier I (see Fig. 1(a)), and kp is the 
Fermi wave length), the effective action S c g is given by 



S a \(j)\ = 



Tui 



2 \ dr 
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j dr j dr'a(T — r') cos 
Jo Jo 



S a [<t>], 



(2) 



In this equation, (3 = 1/ksT, M — C (h/2e) 2 is the mass 
(C is the capacitance of the junction) and the potential 
£/(</>) can be described by 



= - 



dX </>Ij(A</>) - <\> 



(3) 



where Ij is the Josephson current and I ex t is the external 
bias current, respectively. The dissipation kernel a(r) is 
related to the quasiparticle current I qp under constant 
bias voltage V by 



(4) 



at zero temperaturei 24 ' 25 i 26 

Below, in order to investigate the effect of the nodal- 
quasiparticles and the ZESs on MQT, we derive the ef- 
fective action for two types of s/d x junctions, i.e., \ = 
and 7t/4, with \ being the mismatch angle between the 
normal of the insulating barrier (I) and the crystalline 
axis of the d-wave superconductor (see Fig. 1(a)). In 
the case of s/do junctions, ZESs are completely absent. 8 
On the other hand, in s/d^/4 junctions, ZES are formed 
near the interface between I and the d-w&ve superconduc- 
tor d v /4- Note that the influence of the ZES is maximized 
at x — 7r /4 in the range < \ < ir/2. s 

First, we will calculate the potential energy U which 
can be described by the Josephson current flowing 
through the junction (Eq. (3)). In order to obtain the 
Josephson current, we solve the Bogoliubov-de Genncs 
equation with appropriate boundary conditions Then, 



in the case of low temperatures (/? 1 <C Aq) and high 
barriers (zq 3> 1), we get the analytical expressions as 
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where 
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e = A d /A s , K is the complete elliptic integral of the 
second kind, Iq — irA s /2eRN is the Josephson critical 
current for the s/s junctions, Rpj is the normal state re- 
sistance of the junction, Rq = h/4e 2 is the resistance 
quantum, and N c is the channel number at the Fermi 
energy. In this calculation, we have assumed that the 
amplitude of the pair potential is given by A s for s, 
A d cos 26» ee A rfD (0) for d , and A d sin26» ee A rf)r/4 (0) for 
d-w/4, where cosO = k x /kp. By using the analytical ex- 
pression of the Josephson current (5), we get 
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where r\ = I ext /I, 



ext/-«Ci(C2)- 



As in the case of s/s 
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and d/ d 10 i 12 junctions, U can be expressed as a tilted 
washboard potential (see Fig. 1(b)). 

Next we will calculate the dissipation kernel a(r). In 
the high barrier limit, the quasiparticle current is given 
in terms of the convolution of the quasiparticle surface 
density of states (DOS)^ i.e., 



= xEM 5 



dEN L {E,0)N R {E + eV,9) 



x [f(E)-f(E + eV)], (9) 

where ti w cos9/zo is the transmission amplitude of the 
barrier I, N L ^(E,9) is the quasiparticle surface DOS 
(L = s and R — do or GL/4), and f(E) is the Fermi- 
Dirac distribution function. The surface DOS for s-wave 
superconductors is 



N S (E,9) =N S (E) =Re 



\E\ 



(10) 



In the case of do, no ZES are formed. Therefore the angle 
9 dependence of the DOS is the same as the bulk, i.e., 



N do (E,9)=Re 



\E\ 



On the other hand, the DOS for d^/4 is given by 2 ^ 
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The delta function peak at E = corresponds to the 
ZES. Because of the quasiparticle bound state at E = 0, 
the quasiparticle current for the s/d^/4 junctions is dras- 
tically different from that for the s/dg junctions in which 
no ZES are formed^ In the limit of low temperatures 
<C A s ), we can obtain the analytical expression of 
the dissipation kernel a(r) as 



a(r) 
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where K\ is the modified Bessel function. For |r| ^> h/ A s 
the dissipation kernel decays exponentially as a function 
of the imaginary time r, i.e., 
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The typical dynamical time scale of the macroscopic vari- 
able <f> is of the order of the inverse Josephson plasma 
frequency uj p which is much smaller than A s . Thus the 
phase varies slowly with the time scale given by h/A s , 
then we can expand 4>(t) — 4>(t') in Eq. © about r = t'. 
This gives 



h d(j){T) 



2e dr 



(15) 



Hence, the dissipation action S a acts as a kinetic term so 
that the effect of the quasiparticles results in an increase 
of the capacitance, C — > C + SC = C ren . This indi- 
cates that the quasiparticle dissipation in s/d junctions is 
qualitatively weaker than that in in-plane d/d junctions 
in which the super-Ohmic (a(r) ~ |r| 
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or Ohmic dissipation (a(r) ~ r ) 
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appears. At 



zero temperature, the capacitance increment 5C can be 
calculated as 



° 2 Rt 



6C=l 8a%A.£ 



N 
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In this equation 
A(e) = 



X 



x E(y/l - e 2 x 2 ) - 2e 2 x 2 K(Vl - e 2 x 2 ) ,(17) 

where E is the complete elliptic integrals of the first 
kind. Fig. 2 shows the dependence of the capacitance 
increase SC on A s , where we have used Rn = 3.68il 21 




FIG. 2: (Color online) The s-wave gap A s dependence of the 
capacitance increment SC due to the low-energy quasiparti- 
cle dissipation for s/do (dashed line) and s/d n /4 (solid line) 
junctions. 



and Ad — Aybco = 20.0meV. SC is rapidly decreasing 
with increasing A s . Thus, in the case of large A s , the 
influence of the low-energy quasiparticle dissipation on 
the macroscopic quantum dynamics becomes very weak. 
This is a clear indication of quasiparticle-tunneling block- 
ade effects in isotropic s-wave superconductors. As an 
example, if s —Nb and d = YBCO (A s = 1.55mcV, and 
£- 1 = A s /A d = 0.0775), we get SC = 9.5fF for s/d 
and SC = llpF for s/d v n ■ Therefore, in this case, 
the effect of the ZESs on the macroscopic quantum dy- 
namics is considerably stronger than that of the nodal- 
quasiparticles. 

Next, we will investigate MQT in s/d junctions The 
MQT escape rate from the metastable potential (Fig. 
1(b)) at zero temperature is given by^ 3 - 

1 - ■ - (18) 



r 
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By using the Caldeira and Leggett theory^ the MQT 
rate is approximated as 

Wp(fy) 
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is the Josephson plasma frequency (i = 1 for do, i = 2 



for d 



tt/4 



<zi = 1, a-i — 2, and M r , 



(ft/2e) 2 C* ron ) and 



B(jf) = S e s[4>B]/fi is the bounce exponent, that is the 
value of the the action S c g evaluated along the bounce 
trajectory <^b(t). The analytic expression for the bounce 
exponent is given by 



B(rj) 
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where b\ = 12/5 and 62 = 3V2/5. In actual MQT exper- 
iments, the switching current distribution P(j]) is mea- 
sured. P(rj) is related to the MQT rate T(rj) as 



P(r/) = — T(rj) exp 
v 



(22) 



where v = \dr)/dt\ is the sweep rate of the external bias 
current. At high temperatures, the thermally activated 
(TA) decay dominates the escape process. Then the es- 
cape rate is given by the Kramers formula^ 3 - 



2tt 



exp 



Up \ 
k B T 



(23) 



where Uo is the barrier height. Below the crossover tem- 
perature T*, the escape process is dominat ed by MQ T. 
In the absence of dissipative effects, T* cx y/ jc /{C /A), 35 
where jc = I c /A is the Josephson critical current den- 
sity (A is the junction area) and C/A — ej/dj (ej is the 
permittivity of I and di is the barrier thickness). Impor- 
tantly, T* is reduced in the presence of dissipation ! 33 ! 34 
In order to see explicitly the effect of the quasiparticlc 
dissipation on MQT, we numerically estimate T*. We 
determine T* from the relation 



TA( T *\ _ MQT 



a lA (T*) 



(24) 



where a TA (T) and cr M Q T are the standard deviation 
of P(j]) for the temperature-dependent TA process and 
the temperature-independent MQT process, respectively. 
Presently, no experimental data are available for ideal 
highly under-damped s/d junctions with large McCum- 
ber parameters (3m = (2e/h)I c CR^ g » 1 (R sg is the 
subgap resistance). Thus, we estimate T* by using the 
parameters for an actual Nb/Au/YBCO junction^ 1 - ( 
C = 0.60pF, I c = 95.2/xA at % = 0, and R N = 3.680) 
in which (3m ~ 1-5. We also assume Id — Ic 2 — 95.2/iA 
for simplicity, A s = A Nb = 1.55meV, A^ = A Y bco = 
20.0meV^ and vI Cz = 0.0424A/s. In the case of s/d a 
junctions, we obtain T* = 336mK for the dissipationless 
case (C ren = C) and T* = 333mK for the dissipative 
case (CVen = C + SC). Thus, the influence of the nodal- 
quasiparticle is negligibly small. On the other hand, in 
the case of s/d^/4 junctions, we get T* = 601mK for the 
dissipationless case and T* = lOlmK for the dissipative 
case. Therefore, as expected, the ZES have a larger in- 
fluence on MQT than the nodal-quasiparticles. However, 
the T* suppression is small enough to allow experimental 
observations of MQT. 



B. Macroscopic quantum tunneling in 
Nb/Au/YBCO junctions 

In the above calculation, we have assumed that the 
junction is ideal, i.e., the insulating barrier (I) is perfect 
so that we have ignored dissipation by a shunt resistance 



which is caused by imperfections of I. However, in prac- 
tice, it is currently very difficult to fabricate such a per- 
fect junction. Here we will investigate MQT in realistic 
junctions to check the feasibility of the MQT observation. 

Recently, Nb/Au/YBCO ramp-edge junctions have 
been fabricated using the pulsed laser deposition tech- 
niquej 21 ! 37 These junctions can be regarded as hybrid s/d 
junctions between the proximity-effect induced s-wave 
superconductor in Au and the <i-wave superconductor 
in YBCO^i Importantly, the mismatch angle x °f the 
YBCO crystal can be artificially varied as a single pa- 
rameter. By use of these junctions, angle-resolved elec- 
tron tunneling experiments— and scanning- SQUID mi- 
croscopy measurements^ have been performed in order 
to test the momentum dependence of the Cooper pairing 
amplitude and phases in YBCO. 

The main differences between these junctions and the 
junctions with the clean interface discussed in Sec. II 
A are as follows. Firstly, in actual Nb/Au/YBCO junc- 
tions, (i) the ohmic dissipation in the shunt resistance 
will dominate over the intrinsic quasiparticle dissipation 
resulting from the nodal-quasiparticles and the ZESs, be- 
cause the observed McCumber parameter (3m is the order 
of unity. This observation can not be explained only by 
taking into account an intrinsic quasiparticle dissipation. 
Thus we have to consider an extrinsic Ohmic dissipation 
source. Moreover, due to the surface roughnes a 8 ' 36 and 
the very low transparency 3 ^ of the actual insulating bar- 
rier I, (ii) the contribution of the ZESs to the Josephson 
current is expected to be small, so the sin 2<f> component 
of the Josephson current can be neglected. From (ii), the 
actual x dependence of the Josephson current is given by 



Ij{x,4>) = Mx)sin0, 



(25) 



with Ic(x ~ 7r /4) ~ 0£L This can be nicely fitted by 
a theoretical model (Eq. (1) in Ref. [2l|) in which the 
effect of the ZESs is neglected. 

Then the effective action which describes the 
Nb/Au/YBCO junction can be expressed as 



ScS [4 
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2e 
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(27) 



where rj = I ext /Ic(x)- The kernel describing the Ohmic 
dissipation in the shunt resistance can be modeled as^ 4 - 



h Rq 1 
u s {t) = —— j ' 



(28) 



where the shunt resistance R s typically corresponds to 
R S g- Along the same lines of the method used above, one 
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FIG. 3: (Color online) The mismatch angle \ dependence of 
the crossover temperature T* for an actual Nb/Au/YBCO 
junction with (closed circles) and without (open squares) the 
Ohmic dissipation in a shunt resistance R a (inset). 



finds that the MQT rate for the Nb/Au/YBCO junction 
is given by 



2tt 



x exp 



- J B(x,r?) + j— — (1-r? ) 



,(29) 



where the plasma frequency and the bounce exponent is 
respectively given by 



(30) 



B( X ,V) = ^jIc(x)M(l- V 2 )i, (31) 

and £(3) is the Riemann zeta function. The second term 
in the exponent results from Ohmic dissipation in the 
shunt resistance. 

From the MQT rate formula [Eq. ((29)) ]. we numer- 
ically estimate the \ dependence of T* . We have ex- 
tracted Ic(x) horn the theoretical fittin g cu rve of the 
experimental data (see Fig. 3(d) in Ref. [21() and used 
C = 0.60 pF, R N = 3.680, and vl c = 0.0424A/s. In 
this calculation we have assumed R s = Rn for simplic- 
ity. Fig. 3 shows the dependence of T* on the mismatch 
angle x- By changing we can control the macroscopic 
quantum dynamics, i.e., the MQT rate of the junction 
artificially. Note that the maximum T* is attained when 
X = 0. Although T* is reduced by the dissipative ef- 
fects, its magnitude (T* 0.12K at x = 0) is compa- 
rable to high quality s/s junctions (T* w 0.3K)5^ and 



much larger than YBCO junctions (T* w 0.05K)J. This 
can be attributed to the high Josephson critical current 
density jc ~ 14kA/cm 2 and the weak Ohmic dissipative 
effect despite the small (3m- Therefore, MQT should still 
be experimentally observable using the current junction 
fabrication and measurement technology. Moreover, if 
we increase jc or decrease the capacitance per area C/A 
of the junction, we can get still larger T*. 



III. SUMMARY 

In conclusion, MQT in the s/d hybrid Josephson junc- 
tions with the perfect insulating barrier has been theoret- 
ically investigated using the path integral method. The 
effect of the low energy quasiparticles on MQT is found 
to be weak. This can be attributed to the quasiparticle- 
tunneling blockade effect in the s-wave superconductor. 
We also investigated MQT in a realistic s/d junction and 
showed that the expected T* is relatively high in spite of 
the small (3m- These results strongly indicate the high 
promise of s/d hybrid junctions for quantum computer 
applications. 

Finally we would like to comment the advantage of 
the s/d junctions over the s/s and d/d junctions. As 
compared to s/s junction, s/d junctions with the clean 
interface will have much larger I c . So we can expect that 
such s/d junctions show higher T* than the s/s junctions. 
In do/ do junctions, the nodal-quasiparticles give negligi- 
bly small effect on MQTJ£i±i On the other hand, as was 
shown by Fominov et. al.^ the decoherence time of the 
do /do qubit is not enough for practical quantum compu- 
tation. This indicates that the nodal quasiparticles still 
have large influence on the qubit operation. However, 
due to the quasiparticlc tunneling blockade effect in the 
s/d junctions, the decoherence time of the s/d qubit is 
expected to be much longer than that of the do / do qubits. 
Detailed discussion along this line is an interesting future 
problem. 
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